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We obtain a canonical form of a quadratic Hamiltonian for linear waves in a weakly in- 
homogeneous medium. This is achieved by using the WKB representation of wave packets. 
The canonical form of the Hamiltonian is obtained via the series of canonical Bogolyubov-type 
and near-identical transformations. Various examples of the application illustrating the main 
features of our approach are presented. The knowledge of the Hamiltonian structure for linear 
wave systems provides a basis for developing a theory of weakly nonlinear random waves in 
inhomogeneous media generalizing the theory of homogeneous wave turbulence. 



1 Introduction 

In order to analyze the behavior of a general nonlinear system, the first necessary step is to analyze the 
linear system. The classical examples are finite dimensional Hamiltonian systems with coupled degrees of 
freedom. The dynamical behavior of such systems for small excitations, is determined by the form of the 
quadratic part of their Hamiltonians, which correspond to the linear dynamics. A detailed classification 
of the canonical normal forms for such Hamiltonians was given by Galin and it was summarized in one 
of the appendices of Arnold's book [1]. Another class of examples of nonlinear systems, which behavior 
crucially depends on the linear part, is weakly interacting dispersive waves in continuous media that are 
studied in Wave Turbulence (WT). These systems often have analogs among discrete systems, e.g. a chain 
of coupled oscillators, and correspondingly their quadratic Hamiltonians have analogs among Galin- Arnold 
canonical forms. Examples where WT has important applications are surface gravity waves [2], /3-plane 
turbulence in oceans and atmospheres [3,4], internal waves in the ocean [5], weak magnetohydrodynamic 
(MHD) turbulence [6], Bose-Einstein condensate [7,8] and plasmas [9]. Traditionally, WT theory is applied 
to homogeneous systems. The quadratic interaction Hamiltonian with linear interactions describing such 
homogeneous systems is given by [10] 



H 2 = I A k \a k \ 2 dk + -(B k a k a_ k + B k a k a*_ k )) dk, 



(1) 



where a k is a complex-valued field, the bold face denotes a d dimensional vector in R . The linear canonical 
transformation to the new variables b k that was first used by Bogolyubov in 1958 for the system of fermions, 
brings this Hamiltonian to its normal form 

H 2 = j u k \b k \ 2 dk, (2) 

where ui k is a linear dispersive relationship. The formalism of WT significantly enhanced our understanding 
of spectral energy transfer in ocean, atmosphere, plasma, other system of nonlinear waves [10]. Wave 



turbulence deals with weakly nonlinear waves with quasi-random phases. Using WT, one can derive a 
kinetic equation for the wave spectrum, which evolves due to resonant wave interactions. 

In order for the resonant energy transfers to occur, certain resonance conditions need to be satisfied. In 
particular, for the systems dominated by three wave interactions, such as internal waves in the ocean [5] or 
capillary waves [11], these conditions are given by 

^k = ^ki + ^k 2 > 

k = ki + k 2 . 

(3) 

Similar four-wave resonant conditions should be satisfied for the resonant interactions in the four-wave 
systems. 

Often, WT is not spatially homogeneous and its statistical properties vary in space due to a trapping 
potential, inhomogeneous background density or an inhomogeneous velocity field. Examples of such systems 
are the Bose-Einstein Condensate (BEC) in the presence of a trapping potential [8] and an interaction of 
the long "aged" gravity waves with a swell [12]. In general, the idea is to consider a small amplitude high 
frequency perturbation of a large-scale solution of the dynamic equation (e.g. the condensate). The effects of 
this coordinate dependent background solution can most easily be understood using a wave-packet formalism. 
This formalism was first used to approximate the Schrodinger's wave function by a quasi-monochromatic 
wave by Wentzel [13], Kramers [14], and Brillouin [15]. Their initials give the term WKB approximation. 
The WKB approximation is applicable if the wavepacket wavelength I is much shorter than the characteristic 
wavelength L of a large scale solution 

I 

«- z <i. 

The essence of WKB approach is that the wave numbers, characterizing the wavepacket are the functions of 



coordinates. This is due to the distortion of the wave packets by the media, leading to a spatial wavenumber 
dependence. Such spatial wave number dependence may have a dramatic effect on nonlinear resonant wave 
interactions. Indeed, the resonant conditions (3) may now be satisfied only in a finite part of a domain, or 
for a particular wave packet, only for a finite time. 

The goal of this paper is to use the spatially dependent WKB wave packets to find a canonical form 
for the quadratic Hamiltonian for the inhomogeneous systems. This problem can be considered as an 
extension of certain (oscillatory) members from the Galin- Arnold classification of quadratic forms onto the 
infinite-dimensional and continuous space. The physical motivation for such a formulation is that, since the 
Hamiltonian description is natural for WT in homogeneous systems, it lays down a necessary framework for 
generalization onto the inhomogeneous media. 

To begin, we write down the general Hamiltonian for the system of linear waves propagating in the 
inhomogeneous background. We will show in Section 2 that such general Hamiltonian for the variable a q is 
given by the following quadratic form: 



H = j yA(q,q 1 )a cl a^ i + -^B(q,q 1 )aqa-q 1 +B*(q,q 1 )a* q atq 1 jjdqdci 1 . 



(4) 



The main result of the present paper is that this general Hamiltonian (4) can be transformed to the following 
canonical form 

H = / c kx [w kx - x • V x w kx + i{w kx , -}]c kx <ik(ix. (5) 

Here cj kx and c kx are the new position-dependent dispersion relationship and normal field- variable, corre- 
spondingly, and a Poisson bracket is defined by 

{/,<7} = V k /-V x </-V k <?-V x /. 

This novel canonical Hamiltonian is a generalization of the Hamiltonian (2) for the inhomogeneous sys- 
tems. The approach used in the paper can be viewed as a generalization of the Bogolyubov transformation, 



which diagonalizes Hamiltonian (1) to the form (2) via canonical transformation. Similarly, Hamiltonian (4) 
can be transformed into the canonical form (5), however, now using near-canonical transformations. In the 
Hamiltonian (5), the second and the third terms in the brackets correspond to the higher order corrections 
to the dispersion relation due to the inhomogeneity. We prove that the Hamiltonian (4) can be transformed 
into a canonical form (5) in the case when the heterogeneity is weak. Formally, the requirement of weak 
inhomogeneity means that the coefficients A(q, qjj and -B(q, qi) are strongly peaked at q — qx = 0, i.e. 
A(q, q x ) = for |q — q x | > e for some small e. Based on this requirement, we derive below the re-normalized 
dispersion relationship and the transformation formulas from a k to c kx accurate up to the first order in e. It 
turns out that just Bogolyubov's transformation is not enough in this case. The phase coordinate systems 
should also be perturbed by a near-identity transformation in addition to the Bogolyubov's rotation. Then, 
the Hamiltonian becomes diagonal up to the first order in e. 

From the novel canonical form of the Hamiltonian given by Eq. (5), the traditional radiative action 
balance equation can easily be obtained: 

<9rt kx 

— -— + Vk^kxVxWkx - V x u; kx V k n kx = 0, 

or, shorter, 

—K7- + Vkx, ™kx} = 0, (6) 

where n kx denotes the ensemble average of the squared amplitude of the wave, i.e., n kx = (|c kx | 2 ). Equa- 
tion (6) is now a standard equation which is used in statistical modeling of wave systems [16]. It is a wave 
analog of the Liouville's theorem, or the continuity equations for distribution function of statistical me- 
chanics [17, 18]. In wave systems there are wave quasi-particles instead of particles, following different rays 
instead of individual particle trajectories. Wave action distribution function moves in the multidimensional 
wave-number-coordinate phase space. 



The paper is organized as follows. In Section 2, we give simple and instructive examples that motivate the 
study of inhomogeneous WKB systems. In Section 3, we introduce the window transforms and other formulas 
that will be extensively used later. In Section 4, we discuss the case of a nearly-diagonal Hamiltonian. We 
show how it can be transformed to a canonical form (2) and provide a couple of representative examples. 
In Section 5, we study the Hamiltonian in a general form (4). We present the series of near-canonical and 
near-identical transformations that bring the Hamiltonian (4) to the form (5). We also demonstrate the 
application of this approach to the nonlinear Schrodinger equation with the condensate. 

2 Motivation 

In this Section, we motivate our study of inhomogeneous systems by considering an example of the interaction 
of short waves on the background of a long wave. We describe the cases of both three-wave and four-wave 
Hamiltonian systems. 

Three-wave case. The quadratic Hamiltonian of a three-wave system with small scale perturbations 
on the background of the large scale excitations is given by Eq. (4). In order to show that, we start with a 
standard three- wave Hamiltonian [10]: 

H 3 = ft k |a k | 2 dk +- (v£Xaia m + c.c.) ^dktfldm, (7) 

where V is an interaction coefficient. Then, the equations of motion for the variable a^ assume a standard 
form 

™k = — ^ = ^k«k + / (^lmaiOm^im + ( V Li) ai<C<5 km J (Mm. (8) 



Suppose that a large-scale solution of (8) is given by C k . We consider a perturbed solution a k = C k + Ck 
where Ck is a small-scale perturbation of Ck- Equation of motion for Ck attains the following form: 

ic k = O k C k + / -jV\Zn{C\C m + C m Cl + ClCm)^ + (F km J {C\c* m + C^ci + cic^)5 km 

Now, we use the fact that Ck is a known exact solution for Eq. (8) to obtain 



dldm. 



ic k = / A(k, l) Cl dl + - / B(k, ly^dl + / -V^dCmS^ + (vilm)* Qc^, 






dldm, (9) 



where 

A(k,l) = ^ 1 k + y 1 k k _ 1 Ck-i+(y k ,i-k)^i-k, (10) 

s(k,i) = 2y k k :Jc k _i. 

Equation (9) corresponds to the following Hamiltonian 

H = ( A(k, l)cic k (ildk +]- f [B(k, l)cI lCk + c.c] dkdl + - I Vi^cic™^ + c.c.jdkdldm. (11) 

This appears to be a standard form of the Hamiltonian for the wave system dominated by three wave 
interactions in the inhomogeneous media. Quadratic in c k part of this Hamiltonian is given by (4), while 
cubic part of this Hamiltonian is a standard three-wave interaction Hamiltonian. 

Four- wave case. Four- wave systems are similar to three- wave systems when small scale perturbations 
on the background of the large scale excitations are considered. Indeed, we show below that the quadratic 
part of a four- wave Hamiltonian of small scale perturbation has the form (4). We start from a standard 
four- wave Hamiltonian [10]: 

H 4 = / S7 k |a k | 2 dk + - / T^ s alala m a s 5^ s dkdldmds, (12) 

where T is an interaction coefficient. The corresponding equation of motion takes the form 

id k = r2 k a k + / T tns ala m a s 5 ms dldmds. (13) 



We then consider a perturbed solution a k = C k + Ck where c k is a small-scale perturbation. Assuming that 
Ck is an exact solution to the equation of motion with Hamiltonian (13), we obtain the following equation 
of motion for c k 

ic k = tt k c k + / T^ 1 s (2C 1 *C m c s + clC m C s )5™ s dldmds + 

/ T™ s (C?c m Cs + 24c m C s )5™ s dldmds + / T^ s 4c m c s 5^ s d\dmds. 

Since Ck is a known large scale solution, we obtain 



ick = / A(k, s)c s ds-\ — / B(k, s)c*_ s ds 



+ 



2^1mClCm + (W km )* Ci C* m 



dldm 



+ / T™ s 4c m c s 5™ s dldmds, (14) 

where we defined the kernels A(k, s) and i?(k, s) as 

A(k,s) = Q s 5* + 2 i T^CfCmS^dldm (15) 

S(k,s) = 2 J T^ s C m C\5^- s dldm, 

and 

wL = 2 I rJSc:^ ds. 



lm ^ / 1 ml L/ s t 'lm t 

The linear part of equation (14) has the same form as the linear part of the corresponding equation obtained 
for the three- wave case (9). Thus, this linear part corresponds to the same first two terms as in Hamiltonian 
(11). Note also similarity of the quadratic terms in (9) and (14). Note that (14) correspond to the following 



Hamiltonian 



H = J A(k, l)dc£dldk + - f [B(k, 1)^4 + c.c] dkdl 
+ 9 / [^lm c k c l c m + c.c.]dkdldm + 



+ o / T^clctcmCsS^dkdldmds. 



(16) 

This appears to be a standard Hamiltonian for the wave system with four-wave interactions in the presence 
of spatial inhomogeneity. Indeed, the quadratic part (first line) of this Hamiltonian is the Hamiltonian (4). 
Cubic term is the three-wave interactions with the background large scale wave (i.e. four wave interaction 
where the role of the fourth wave is assumed by the background wave) . Notice that unlike traditional three 
wave interactions in a homogeneous environment, momentum is not conserved by this term. This is the 
effect of breaking of spatial symmetry by an inhomogeneous background. Lastly, the quartic term (third 
line) is the standard four wave interactions Hamiltonian. We show in this paper that the quadratic part of 
this Hamiltonian may be reduced to the novel canonical Hamiltonian for spatially inhomogeneous systems 
(5). 

In this section we have demonstrated that if the general wave system is dominated by three-wave or 
four-wave interactions, and consists of short scale waves superimposed on known large-scale motion, its 
quadratic Hamiltonian is given by the Eq.(4). 

3 Preliminaries 

In this Section, we set up the stage for formulation of our results. Here, we give basic definitions, and obtain 
frequently used formulas. 



We use the following definition of direct and inverse Fourier transforms: 



<7(x) = Jg(ky k *dk. 



Next, we generalize the Fourier transform to spatially inhomogeneous systems. In order to do that, we use 
a window transform of g(x): 

r[5(x)] =5(x,k) = ^i^y/( e *|x-xo|Mx ) e -^ Xo dx . (17) 

Here, f(x) is an arbitrary fast decaying at infinity window function. The parameter e* is defined by the 
spatial scales of the inhomogeneity and the propagating wave-packets in the following manner. First, we 
introduce the characteristic length of inhomogeneity to be of the order 1/e. Then, we take the width of the 
window, which is of the order 1/e*, to be much smaller than the characteristic length of inhomogeneity. On 
the other hand, the width of the window is chosen to be much larger than the wavelength of the waves that 
propagate in the inhomogeneous medium, which is of the order 1. Therefore, we have 

e < e* < 1. (18) 

The special case when f(x) = exp(— x 2 ) is called Gabor transform [19]. Note that, when e* approaches zero, 
f(e*x) approaches the constant function with the value one. Consequently, the Gabor transform becomes 
a Fourier transform. Therefore the Fourier transform can be seen as an averaging over an infinitely large 
window. 

The inverse of the window transform (17) is given by 

g(x) = J g^e** dk, (19) 

where we have used /(0) = 1. We emphasize that Eq. (19) and all the formulas that we obtain below can 

10 



be obtained using any fast decaying at infinity window function and are independent of the particular form 
of fix) as long as it is sufficiently smooth. 

Now, we present the formulas for the window transform, which will be useful later. First, we express 
the window transform g(x, k) in terms of the Fourier transform g(k) 

fif(x, k) = JL | /(|k - q|/e*)e^- k )- x 5(qMq. (20) 

Next, we express the Fourier image g(k) in terms of the window variable g(x, k) 

g(k)=(^=] |s(x,k)dx. (21) 

By combining Eqs. (20) and (21), we obtain the following formula 

g(x,k)=(±=} |/(|q-k|/ £ *)e^- k )' x 5(x / ,q)dqdx'. (22) 

After introducing notations and formulas that will be extensively used below, we proceed to the discussion 
of the main results of the paper. 

4 The case of nearly- diagonal Hamiltonians. 

4.1 Formulation and Proof of the Lemma 

Let us start with Hamiltonian (4) without off-diagonal terms, so that -B(q, qi) = 0. This is a typical 
Hamiltonian for linear waves in weakly inhomogeneous media [10] expressed in terms of Fourier amplitudes 



a q and a* as 



H= / fi(q,qi)o q o^ 1 dqdq 1 , (23) 



11 



with a Hermitian kernel 0(q 1; q) = fi*(q, qi), which is strongly peaked at q — qj = (0 has a finite support 
around q ~ q x ). Therefore, we subsequently assume that there is a small parameter e for which 

fi( qi - q) = 0, (24) 

when |q — qi| > s. A particular choice 

^(qi,q) =w(qi)5(q-q 1 ), 

leads to the familiar form of the Hamiltonian (2). 
The equation of motion for a^ is 

r)ff I' 
»Ok = "Fir = / ^qkaqdq- (25) 

Lemma. Consider the Hamiltonian (23) with £l(q, q\) being a peaked function of (q— q{) (satisfying (24)) 
and a smooth function of {q+ q 1 ). Then there exist a near- canonical change of variables a k — ► a kx such that 
in the new variables the equation of motion can be written in the Hamiltonian form 



.d_„ 8H 



.f 



'dt^~ 8a* k J 

with the filtered Hamiltonian in a canonical form 

H f = J a kx [. kx - X • V x . kx + i {wtaj .}]alJkd X , (26) 

where uj kx is the position dependent frequency related to Q(q, qi) via the Wigner transform 

u k x = I e imx n(k - m/2, k + m/2) dm. (27) 

Proof. In order to obtain the new variables a^x) we first make a window transform via (17), which is then 
followed by a near-identical transformation. The idea of the proof is to use the peakness of the kernel 
H(q 1 ,q). We make a Taylor expansion around the peak and then by neglecting the higher order terms we 

12 



obtain the desired result. 



We make a window transform from a^ to a^x using Eq. (20). Differentiating Eq. (20) with respect to 
time, using the Eq. (25), and applying the inverse formula (21) yield 

4«kx = (p) y/dk-qil/e^e^-^-^q^Oqdqdq! 
Let us change variables from (q, q x ) to (p, m) as 



(28) 



q = p-m/2, (29) 

q 1 = p + m/2. (30) 



Below it will be convenient to use 



F(p,m) = fi(p-m/2, p + m/2). (31) 

Next, we will approximate the RHS of Eq. (28) by a variation of some filtered Hamiltonian Hf, i.e., by 
SHf/Sa^ x . We can rewrite Eq. (28) as 



*^kx= (4=) y , /(|k-p-m/2|/e*)e^ +m / 2 - k ) x F(p,m)a p _ mAxi dpdmdx 1 . (32) 



'dt 



Let us make another change of variables p — ► p + m/2 

d 



i^«kx = (-^) y /(|k - p - m| / £ *) e < (p+m_k) ' x ^(P + m/2,m)a p , Xl dpc/mdx 1 . 



(33) 



In order to simplify Eq. (33), we are going to use the fact that ^ qqi and /(k) are peaked functions of (q x — q) 
and k, respectively, and fast decaying at infinity. We also keep only first order terms in spatial derivatives, 
neglecting second and higher order terms. Then, we could write 

/(|k - p - m|/e*) = /(|k - p|/e*) + m • V p /(|k - p|/e*) + h.o.t. (34) 

13 



Similarly, we obtain 

F(p + m/2, m) = F(k + p - k + m/2, m) = F(k, m) + (p - k + m/2) • V k F(k, m) + h.o.t. (35) 

where h.o.t. denotes higher order terms. Now, we substitute the expansions (34) and (35) into Eq. (32), 
and after ignoring higher order terms, we obtain 

d ( i xd 



at \ \iix 



x J (/(|k - p\/e*) + m • V p /(|k - p|/ £ *)) e *(P+— k)- 
x (F(k, m) + (p - k + m/2) • Vki ? (k, m))a P;X1 <ip<im<ixi. 

(36) 

Note that, here we have an expansion with two different small parameters e and e*, which obey Eq. (18). 

In Appendix A, we show how to simplify the RHS of Eq. (36). As a result of this simplification, we 

obtain 

d 
i7r:«kx = ^kxOkx - x • V x w k xakx + *V x u;kx ■ V k akx - «V k ^kx ■ V x a kx + «(V k ■ V x )w kx (x • V x a kx ) 
at 



+i/2(V k • V x )a; kx akx + (Vk ; V x )^kx(Vk ; Vx)akx - (37) 

Here we have underlined the terms that have two spacial derivatives. In the spirit of WKB approximation 
these terms will later be neglected. This equation can be written in the Hamiltonian form 



d 5H 

^TTTOkx 



/ 



at <5a* x ' 

where the filtered Hamiltonian takes form 

H f = / (^(^kx-x- V x w kx )|akx| 2 + (V k - V x )u;k x akx( v k • V x )a kx 
+ ««kx( v xWkx • V k a kx - Vk^kx • V x a kx + 
(V k • V x )w kx (x • V x a kx ) + l/2(Vk • V x )w k xakx)j dkdx. (38) 

14 



Now, we will use the general method of the WKB approximation. We will only keep the terms, which 
are of the first order in a small parameter e. In our case, the small parameter e characterizes the rate of 
spatial change of the position dependent frequency Wkx and the dynamical variable Ok x - To apply the WKB 
approximation to Eq. (37), we neglect the terms that have two derivatives with respect to x (underlined) 
because each spatial derivative is of the order e small. As a result, we obtain the following equation of 
motion 

d 
J7-akx = ^kx^kx - x • Vx^kxflkx + »V x w k x • V k akx - zV k ^kx ■ V x a kx + V 2 ( v k • V x )wk x akx- (39) 

However, Eq. (39) becomes non-Hamiltonian. Indeed, the corresponding functional 



Hf =-- J ((wkx -x- V x u; kx )|a kx | + ia kx (V x w kx • V k a kx - V k u; kx • V x a kx + 

l/2(V k • V x ) Wkx S kx )) dkdx, (40) 



is not self-adjoint if (Vk ■ V x )w kx 7^ 0. Therefore, in order to obtain canonical equations of motion, another 
near-canonical change of variables needs to be performed 

»kx(i) = Skxakx(i), (41) 

where s kx is some time-independent function to be determined below. Note that transformation (41) is 
canonical if and only if 

l^kxl 2 = 1- 

Therefore, we need to find such Sk x that the system becomes Hamiltonian in terms of new variables 6k x and 
the transformation (41) is near-canonical, i.e., |skx| ~ 1- We substitute Eq. (41) into Eq. (39) to obtain 

d 
isu. x oi«kx = Skx [(wkx - x • V x u; kx )a kx + «V x a;kx • V k a kx - ^V k u; kx ■ V x a kx ] 

+a kx [iV x w kx • V k «kx - «V k w k x • V x s kx + i/2(V k ■ V x )u; kx s kx ]. 

15 



If we find Skx that satisfies the equation 

Vk^kx • V x s kx - V x ^kx • V k Skx = 2 S kx(V k • V x )w k x, (42) 

then the equation of motion in the new variables dk x takes the canonical form 

d 

i ol^kx = (^kx - x • V x a; kx )akx - «V k ^kx ■ V x d kx + «V x a;kx ■ V k a kx , (43) 

with the corresponding Hamiltonian (26). In order to find a solution of Eq. (42), we make a change of 
variables 

5kx = 21ns k x, (44) 



to obtain 

Vk^kx • V x gf k x - V x cj kx • V k 5kx = (V k • V x )w k x- (45) 

We find the solution of Eq. (45) using the method of characteristics. The characteristics (x(r),k(r)) are 
given by the following equations 

-^ = V^kx, (46) 

elk 

~T — -v x u;k x , 

CLT 

where r is a parameter along the characteristics. Physically these characteristics correspond to the trajec- 
tories (rays) of WKB wavepackets in the (x, k) space. The solution of Eq. (45) is given by 



<7(x(r),k(r))= f (V k • V x )w kx dr' . 
Jo 



Now, we use Eq. (44) and then Eq. (41) in order to obtain the new variable dk x out of the Gabor variable 
ak x . The dynamics of the new variable dk x is described by the filtered Hamiltonian (26). □ 
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Note that for the special case (Vk ■ V x )cJkx = 0, the Gabor variables cikx provide a Hamiltonian struc- 
ture (26). Then, in this special case, there is no need in performing the second transformation (41). We will 
see below in the examples that this observation may significantly simplify the applications of the Lemma. 

To describe the statistical properties of spectral energy transfer in the systems, it is convenient to define 
a position dependent wave action as 

n kx = (|akx| 2 )- (47) 

Using this definition, one obtains from (26) the familiar form (6) of the kinetic equation (sometimes called 
radiative balance equation) in a weakly inhomogeneous media, by using the Eq. (43) . The resulting equation 
for the time evolution of n^x is the Eq. (6) of the introduction. This is the so-called waveaction transport 
which is typical for WKB systems. 

4.2 Relation to the Wigner Transformation 

One can also derive waveaction transport equation (6) directly from equation of motion (25), without 
obtaining first the Hamiltonian structure (5). To do this we define the Wigner waveaction by using the 
Wigner transformation: 

^ = /e imx (a k+ i m a*_, m )dm. 

(48) 

Then the Wigner waveaction nj^. obeys the same kinetic equation (6). The prove can be found, for example, 
in [20] . The idea of the proof is to calculate the time evolution of waveaction n^ x by using definition (48) 
and equation of motion (25). Then one uses transformation similar to (29,30), and expands n^ x using the 
smallness of m, and finally uses integration by parts to obtain (6) - see [20] for details. 
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To address the question on how waveaction n^ x , defined through the Wigner transformation, is related 
to the wave action nk X) defined through the Gabor variables, we substitute (21) into the definition of (48) 
to write 

2d 

n 



Sc = I 4^ I / (ak +m /2, x ial- m /2,x") eim '* dx ' dx " dm - 



Taking into account that ak+m/2,x' an d ^k-m/2 x" are s ^ ow functions of x' and x", one can neglect this slow 
coordinate dependence relative to fast coordinate dependence in the exponent. This allows to perform x' 
and x" integrations to obtain a delta- function with respect to the m argument. Consequently we obtain 
that these two wave-actions, (48) and (47) are approximately proportional to each other: 

w 

n k,x « n k,x- 

There are several important advantages of our method. First, it allows to rigorously write the equation 
of motion for the field variable ak, x i n addition to the transfer equation of (6). In addition, our approach 
shows how to derive the kinetic equation (6) to a much broader class of nonlinear systems, those described 
by equation (4) with nonzero value of i?(q, q 1 ) 7^ 0, as we show in the next section. Lastly, Hamiltonian 
formulation helps to rigorously establish a wave turbulence theory and to take into account nonlinear wave 
interactions. 

4.3 Example: Linear Schrodinger equation 

Consider a one-dimensional example of a linear Schrodinger equation with a slowly varying potential. This 
equation is also referred to as linearized Gross-Pitaevsky equation. It is used to describe a formation of the 
BEC. It is given by 

iiP = _v^ + U(x)ij). (49) 
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This equation can be written in a Hamiltonian form with the Hamiltonian given by 

H = J (\V X ^\ 2 + U{x)\^\ 2 ) dx. (50) 

In the Fourier space, Eq. (49) becomes 

i^k = k 2 ^ k + f U{k - k^ih) dk u (51) 

with a corresponding Hamiltonian 

H= Q(k,ki)il;i,ipl 1 dkdki, 

where Q,(k,k\) = k 2 b\ + U{k\ — k). Now we can apply Lemma and find that the position dependent 
dispersion becomes 

v kx = k 2 + U(x), 

and the corresponding Hamiltonian in terms of Gabor variables takes the canonical form (5): 

H, = / feh. - *V«« + **., ■>]&*«*■ (52) 

It follows from the Lemma that the Gabor variables provide a canonical description of the system (49). 
Indeed, since VfcVx^fcx = 0, we do not need to make a near-identity transformation (41) in the Lemma. 
Therefore, it is instructive to obtain the same result by directly applying Gabor transform to the both sides 
of Eq. (49). We have 

T[Vlip] = Vl^kx + 2ikV x tjj kx ~ k 2 i>kx ~ 2ikV x i) kx - k 2 tjj k x- (53) 

To obtain this equation we have neglected the higher order derivative of the Gabor variable, since it is slowly 
varying in x. Next, we use the linear expansion of the potential U(xq) pa U(x) + (xo — x)V x U{x) to find 

T[U(x)1>] = (U(x) - xV x U(x))^ kx + iV x U(x)V k 4> kx . (54) 
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Combining Eq. (53) with Eq. (54), we obtain 

d ~ ~ 

i-^2^kx = (k 2 + U(x) - xV x U(x))tp kx + iV x U(x)V k ip kx - 2ikV x tp kx . 
at 

Hamiltonian that corresponds to this equation is given by Eq. (52). 

4.4 Example: an advection-type system 

Let us consider an advection-type system. For simplicity of calculations let us restrict our attention to a one 
dimensional case, although a general dimensional case can also be considered. An advection-type system 
has a Hamiltonian of the form 

H = i J U{x) [ip(x)V x ip*(x) - i/>*(x)V x i/>{x)] dx. (55) 

with the corresponding equation of motion 

i^{x) = -iV x (U(x)ip(x)) - iU{x)V x il>{x) 

= -i(2U(x)V^(x)+V x U(x^(x)). (56) 

In the Fourier space, this system is described by the Hamiltonian 

H= Q(k,ki)'ii) k 'ijj k ' 1 dkdq kl , 

with the kernel 

n(k,k 1 ) = (k + k 1 )U(k 1 -k), (57) 

After applying Lemma to Hamiltonian (55), we obtain the following canonical form 

b, = / JUu*. - xv«« + i{o, ta , m^ x . (58) 
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where ip are the new variables. 

oj kx = 2kU(x), (59) 

is a position dependent frequency. Note that in this case we have V^V^Wfa 7^ and the near-canonical 
change of variables given by Eq. (41) had to be performed. 

We can also obtain the same result by directly applying the Gabor transform to Eq. (56). Using the 
slow dependence of U(x) on x (disregarding the second derivative and higher), we obtain 

T[U{x)V x i>{x)\ « 
(U(x) - xV x U{x) + iV x U(x)V k ) {V x + ik)i> kx . (60) 

Similarly, we have 

T[V x U(x)^(x)] « V x U{x)i> kx . (61) 

Substituting Eqs. (60) and (61) into Eq. (56), we obtain 



2iV x U(x) + 2iV x U(x)kV k - iV x U{x))^ kx . 



(62) 



Using Eq. (59), we rewrite Eq. (62) as 



i-Kii'kx = {u kx - xV x iv kx + V k V x uJ kx V k V x + i(V x uj kx X7 k - V k uj kx V x ) + 

ixV k X7 x uj kx X7 x + -iV k V x uj J ip kx . (63) 

As in the Lemma, we neglect the higher order terms with the two derivatives over x (underlined in Eq. (63)). 
In order to obtain the canonical form of the equation of motion, we need to make a near-canonical trans- 
formation 

i>kx = ftpkx, 
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where / satisfies 



^k^kx^xf - VajWfcccVfe/ = -/VfcV x W fcx - (64) 

For this special case, we obtain 

U(x)V x f = V x U(x)(^f + kV k fY (65) 

We have to find the solution for Eq. (65) such that / — ► 1 when X7 x U(x) — ► 0. Therefore, we need to find 
a solution in the form / = 1 + g, where g satisfies \g(k,x)\ <C 1. Let us try to find a solution in the form 
/ = f(x), i.e., independent of k. Then we have 



/(x) = Cy/Uix), 
where C is an arbitrary constant. We expand U(x) around some point of reference xo as 

U(x) « U(x ) + (x- x )V x U(x ). 
Let us choose the constant to be C = l/y / f/(xo) then we obtain 

/wl+ 2 (a: - aro) -^r- 

If Va;C/(xo) ~ s and |x — xq\ <€. 1/e then / ~ 1 and the transformation is near-canonical. 

5 General case of waves in weakly- inhomogeneous media 
5.1 Theorem 

The result of the Lemma can be generalized onto a much broader class of Hamiltonian given by (4). Thus 
let us consider (4) with both A and B being peaked functions of q — q l5 which corresponds to waves on a 
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weakly inhomogeneous background. Note that 

A(q 1 ,q) = A*(q,q 1 ), (66) 

because the Hamiltonian is Hermitian. Moreover 

B(-q,-qi)=B(qi,q). (67) 

Condition (67) does not really restrict our choice of the coefficient S(q, q x ). Indeed, we can consider any 
function B and then represent it as a sum of two components 

B(q,q 1 ) = B'(q,q 1 )+B"(q,q 1 ), (68) 

where 

B'(q, qi ) = i(B(q, qi ) +B(-q 1 ,-q)), 
S"(q,qi) = ^(B(q,qi) - B(-q l5 -q)). 

When we plug Eq. (68) into Hamiltonian (4), the part of the integral with B vanishes. 

From now on, we will omit indices wherever it does not lead to confusion and denote a = d^ and 
a_ = a_k,x- Also, we introduce some convenient notations, which we will use throughout the rest of the 
paper. For any function <£>(k), we denote its even part as (p ev and its odd part ip 0( i 

Vev = -{ip + f-), 
"Pod = ^((f-(p-). 

We now are ready to formulate the main theorem of this paper: 

Theorem. Consider a Hamiltonian (4) with A(q, q 1 ) and B(q, g x ) being peaked functions of (q — q\) with 
the same parameter e, i.e., A(q, q ± ) = and B(q, q{) = when \q— q ± \ > e. Suppose that conditions (66) 
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and (67) are satisfied. Let us introduce new notations 

H = f A(k-m/2,k+m/2)e imx dm, 

A = f B(k-m/2,k+m/2)e imx dm, 

v = »[A], 

v = 9f[A]. 

Suppose that Hamiltonian (4) has a dominant diagonal part, i.e. fj, ev > v and v = O(e). Then, there exists 
a new- canonical change of variables from ak to Ck x o,nd the evolution of system (4) can be approximately 
described by the following equation of motion 

d _ 6H f 

with a filtered Hamiltonian 

Hf = / Ckx[u — x ■ V x u> + i{uj, -}}c kx dkdx. (69) 

The position dependent frequency is given by the formula 



The transformation from the Fourier variables a k to new variables c kx is given in the proof. 
Proof. The proof consists of three main steps. In order to diagonalize Hamiltonian (4), we 

1. apply the Lemma to simplify the Hamiltonian using the peakness of the kernels, 

2. perform Bogolyubov transformation to diagonalize the O(l) part of the Hamiltonian, 

3. make a near-identity canonical transformation to diagonalize the 0{e) part of the Hamiltonian. 
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Step 1: Applying Lemma. 

Similarly to Eq. (26), we can write a filtered Hamiltonian for Eq. (4) as 

Hr = / a(fi — x • V x/ u + i{fi, -})a*<ikdx + 

- / [a [A — x • V X A + i{\, -}]a_dkdx + c.c], 

(70) 

here, as usual c.c. stands for complex conjugate. From property (67) and definition (69) it follows that 

A(-k,x) = A(k,x). 

Step 2: Bogolyubov transformation. 

In this step we apply the usual Bogolyubov transformation. Before doing that notice that Hamiltonian (70) 

consists of two parts 

o-(i) _ o-(i) . o-(i) 
H f ~ H f,l +H f,e' 



where 



Ht^ = / /j,\a\ dkdx. -\ — / i/[aa— + a* a*L]dkdx., 



o-(i) _ o-(i) rj-(l) 
H f,e ~ H f ~ H f,l> 



(71) 



are correspondingly O(l) and 0(e) parts of H\ . In Step 2, we diagonalize the O(l) part using the following 
linear transformation 

a = u kx b + v kx b*_. (72) 

It was shown in [10] that transformation (72) is canonical if the following conditions are satisfied: 

KkxP - l^kxl 2 = 1, 
^kx«-k,x = ^-k.x^kx- 
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Let us follow [10] and choose 



u kx = cosh(£ kx ), 
v kx = sinh(£ kx ), 



(73) 
(74) 



where £ kx is real and even, but otherwise arbitrary function. Then under change of variables given by 



Eq. (72), Hf\ becomes 



H 



(i) 
/■i 



/ n cosh 2 (£) + n_ sinh 2 (£) + 2v sinh(£) cosh(^) \b\ 2 dkdx 



+ 2 



(fi + ft-) sinh(£) cosh(£) + i/(cosh 2 (£) + sinh 2 (£)) (66- + b*b*_)dkdx 



Denote the expression in square brackets multiplying \b\ 2 as to: 



uj = ficosh (£) + \X- sinh (£) + 2^sinh(^) cosh(£). 



Using trigonometric formulas for hyperbolic functions, we obtain 



uj = jUeuC0sh(2£) + fi od + z/sinh(2£). 



r(l) 



In order to diagonalize i?W, we require that the following condition in satisfied 



[i ev sinh(£) cosh(£) + v{ cosh 2 (£) + sinh 2 (£) ) = 0. 



This condition is equivalent to 



tanh(2£) = . 



Since fi ev > u, we can choose cosh(2£) to be positive and, therefore, we have 



cosh(2£) 
sinh(2£) 



Me 



VM 



2 _ v 1 
ev " 



Vi4v - v 2 



(75) 



(76) 



(77) 



(78) 
(79) 



Ki 



In Appendix B, we find the expression for £. Resolving Eq. (75) together with Eqs. (78) and (79), we obtain 

W = LOev + ^odi 
^od = Mod) 



Therefore, we have diagonalized H \ ^ to the form 

Hfl = J bu;b*dkdx. (80) 

Next, we consider the O(e) part of the filtered Hamiltonian. In Appendix C, we show that Bogolyubov 
transformation (72) transforms H}' to the form 

Hfl = I b(-x-V x uj + i{io,-})b* + (abb- + ^b{<p,b-} + c.c^\dkdx, 

(81) 



where 



2v \ \x ev 2 2 



L v 2 uj ev 

V r'ev H'ev 



Combining Eqs. (80) and (81), we finally obtain Hamiltonian in the form 

Hf = [b(u-x-V x u + i{w r })b* + (abb- + ^b{<p,b-} + c.cAdkdx. (82) 

Step 3: Near-identity transformation. 

In Step 2, we diagonalized O(l) part but not all of the 0(e) part. In order to diagonalize complete Hamil- 
tonian, we use the near-identity transformation. This near-identity transformation changes variables from 
frkx to Ckx by the following rule 

6kx = c kx + akcl kjX + A<{7k, c!. kjX }, (83) 
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where, we assume that pk and 7k are O(l) terms and Qk and (/3k{7k, c _k)) are 0(e) which makes our 
transformation indeed near identical. Note that a^, /?k, and 7k are functions of both k and x. Nevertheless, 
for simplicity of notation, we do omit the dependence on x, since it would only unnecessarily pollute the 
notations. In Appendix D, we derive the canonicity conditions for transformation (83). In turns out that 
transformation (83) is canonical if the following conditions are met 

A< = /3-k, 
7k = 7-k, 



a od = -^{jk,Pk}- 



(84) 



Among the coefficients a^, /3k and 7k that satisfy the canonicity conditions we have to choose those that 
will diagonalize the 0(e) part. In Appendix E, we show that such coefficients become 

C k /^k ( Uev 

«k = w«i, 

Pk — 



2uuj ev ' 

7k = • 

(85) 

Note that these conditions are in full correspondence with the canonicity conditions (84). The Hamiltonian 
in new variables up to 0{e) order is 

Hf = / c[uj — x • V x u; + i{u>, -}]c*dkdx. 

This completes the proof of the main result of this paper. □ 
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5.2 Example: Nonlinear Schrodinger Equation — waves on condensate. 

Bose-Einstein Condensate (BEC) is a state of matter that arises in dilute gases with large number of particles 
at very low temperatures [21-26]. BEC can be described by the Nonlinear Schrodinger equation (also known 
as Gross-Pitaevskii equation [27]). Here, we apply the theorem to this well studied model. The evolution of 
the state function ip is described by the following equation 

with a corresponding Hamiltonian 

The term n(t)ip is introduced for convenience as will become clear later. Following [7], let us consider the 
amplitude-phase representation of the order parameter tjj: 

V> = Ae*. 
Now, we introduce Hamiltonian momentum 

p = 2A<p, (86) 

and rewrite Eq. (86) in terms of new canonical variables A and p as 

m - f . m 

5H 

pt = -Ja> 

where 



/ ( (va) 2 + \a* - «(t)A 2 + \ (v P - ^y) d* 



H = / (VA) 2 + ^-A 4 - K (t)A 2 + ± (Vp - mp) dx. (88) 
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Let us consider weak perturbations on background of a strong condensate, 

A = A (o) +A d) j p = p (0) + p(i) j |^(i)|«|^(o)|. (89) 

We now choose n(t) = (A^°>) 2 which gives us p(°> ~ e. Substituting Eq. (89) into Eq. (88) we have 

H = H o + H 2 + H 3 , 

where the subscripts denote the order of the term with respect to perturbation amplitudes. Since in this 
paper we study the linear dynamics, we only consider the quadratic part of the Hamiltonian 



Ho 



J ((VAW) 2 + (A^(A^f + \(V P f + \^V P ■ VA^ + \ P (Vln^°)) • Vp\ Ac. 



Here, we used the fact that the spatial derivative adds one order in e and we neglected the terms of the 
order two and higher. 

In order to apply the theorem to H 2 we first transform to Fourier space and then switch to normal 
variables. Let us denote R = (A^) 2 , S = pW/A<® and T = V X lnA (0 \ We have R = O(l) and S, T = 0(e). 
Transforming H 2 into Fourier space we obtain 

H 2 = f ( (ki • k 2 6l + R 2 -i)A 1 A* 2 + ^ 2 _iki • k 2Pl A* 2 + i(ki ■ k^ 1 - T 2 _iki • (k 2 - ki))pu^ dl2, 

where we used the following simplified notations: 1 = ki, 2 = k 2 and subscript 2 — 1 = k 2 — ki. Next, we 
switch to normal variables using the transformation 

A k = _(a k + al k ), (90) 



1 

V2' 



Pk = 7^(ak-al k ) 



In normal variables, H 2 reads 



H 2 = J I (^k 1 -k 2 5 1 2 + R 2 . 1 + ^T 2 ^ 1 (k 2 -k 1 )Aa 1 a* + 

+ - f-ki ■ k 2 5\ + i? 2 _i - ^S 2 -iki • k 2 + -r 2 _iki • (k 2 - ki) J aia_ 2 + c.c. J dl2. 
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The only part of the coefficient in the second parenthesis we are interested in is the one that satisfies Eq. (67): 

3 i 1 

-ki • k 2 5\ + #2-1 - -S 2 _iki • k 2 + -T 2 _i(k 2 - ki) 2 . 

Since A^°> and p^' are slowly varying functions of x, so are -R(x), S(x) and T(x). Therefore, their Fourier 
transforms are peaked around zero making the terms proportional to T 2 _i(k 2 — ki) 2 of the second order in 
e, which can be neglected. Finally, we can write down the Hamiltonian in the form given in Eq. (4) 



H 2 = / (A{k 1 ,k 2 )a 1 a 2 - + -{B(k 1 ,k 2 )aia* 2 + c.c.)jdl2, 



(91) 



where 



A(ki,k 2 ) = |ki • k 2 5\ + i? 2 _! (92) 

3 i 

B(ki,k 2 ) = -ki • k 2 ^2 + R 2 -i - -S^-iki • k 2 . 

In terms of window transformations, which we denote here as a the Hamiltonian reads 

Hf = / a(fi — x • VfJ, + i{n, -})a*dkdx. + - I [a [A — x • VA + i{\, -}]a-dkdx. + c.c], 

where 

(i = i e im -*A(k-m/2,k + m/2)dm = jk 2 + R(x), (93) 

A = 1 / , e imx S(k-m/2,k + m/2)dm = ^k 2 + J R(x)-^k 2 5(x) + ik-V5(x). 

Up to the first order in e, we have 

v = |k 2 + i?(x), (94) 

v = -\k 2 S(x). 

Here, fi is an even function of k which means that n ev = [i and [i 0( i = 0. Then, the position dependent 
frequency of the small perturbations in the presence of the condensate becomes 



u = vV _ v 1 = \k\y/R(x) + k 2 . 
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Bogolyubov's transformation, a = ub + vb*L, is given by the following coefficients 



H 



5k 2 + 4R(x) 



V7* 



2 _ „2 



l|k| A /k z +i?(x) 
3k 2 + 4i?(x) 



V^ 2 " ^ 2 4|k|^/k 2 + i?(x) 

In terms of variables 6 the Hamiltonian takes the following form 

Hf = / b(u - x • V x w + i{u, -})6*dkcix + ( / abb^dkdx +- 



b ^W,b-} 



dkdx + c.c. , (95) 



where 



a 



<P 



^ „ v 2 i~ k 2 



x-VR(x) 



+ iS{x) 



k\Jk z + R(x) 



2 4|kU/k 2 + J R(x) 



/i- 



5k 2 + 4R(x) 



Finally, we perform the near- identity transformation b = c + ac*_ + /3J7, cl} where 



x- V-R(x) 



o 



/S 



7 



i|k|5(x) 



4(k 2 + i?(x)) 4 A /k 2 T^ 



i(5k 2 + 4 J R(x)) 5 



8|kU/k 2 + J R(x)(3k 2 + 4i?(x)) 



4|k|yk 2 + J R(x) 
5k 2 + 4i2(x) 



The resulting Hamiltonian attains the canonical form (5). 



6 Conclusions 



We have studied the dynamical behavior of the linearized spatially inhomogeneous Hamiltonian wave sys- 
tems. The canonical transformation from the Fourier variables to the new spatially-dependent variables is 
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found for the general class of the quadratic Hamiltonians. In the new variables, the linearized dynamics 
is governed by the canonical diagonal Hamiltonian with the spatially-dependent dispersion relation. The 
waveaction transport equation which coresponds to this Hamiltonian has form (6) which is typical to WKB 
formalism. It was previously obtained for some specific examples, e.g. in plasmas [29] and geophysical 
waves [16, 28]. In this paper, we have given several representative examples illustrating the general results, 
such as the Nonlinear Schrodinger Equation without and with condensate and an advective-type system. 
Further possible areas of application of this formalism include water waves on lakes with variable depth 
or/and presence of variable mean flow, internal waves in media with variable stratification, plasma waves 
on profiles with variable density, geophysical waves in media with variable background rotation rates, etc. 

The new Hamiltonian formalism that is presented in this paper should be crucial for extending the 
WT theory to the spatially inhomogeneous systems. In the spatially homogeneous systems, quadratic term 
in the Hamiltonian corresponds to the first term in Eq. (69). Effect of space inhomogeneity leads to the 
appearance of the derivative terms in the Hamiltonian, which correspond to the slow dynamics along the rays 
in the (k, x)-space. This effect will lead to an interesting interplay of inhomogeneity and non-linearity in wave 
turbulence systems. More specifically, linear dispersion relation becomes spatially dependent. Consequently, 
the resonance conditions change as waves propagate through inhomogeneous environment. As a result of 
this, waves will remain in resonance for a limited amount of time, or the members of the resonant triads will 
change from position to position. The physical implication of this effect may be the weakened flux of energy 
or other conserved quantities through the wavenumber space. Other effect may be an effective broadening of 
the resonances, as resonances will be altered from place to place, so the wavepacket propagating through the 
inhomogeneous environment will be affected by averaged dispersion relation. Another potentially interesting 
effect is an effective three- wave interactions in a four- wave weak turbulence systems, where the role of fourth 
wave is played by inhomogeneity. 
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In order to develop a Wave Turbulence theory for spatially inhomogeneous systems, the kinetic equation 
has to be obtained for the cases with such finite-time wave resonances. This is an exciting task for the future 
work. 

Acknowledgments Y. L. was supported by NSF CAREER grant, DMS 0134955. 

A Window transform 

Let us consider the RHS of Eq. (33) term by term. 

1) foF 

(-?-) J /((k - p)/e*)e^ +m - k ^F(k, m)a p ^dpdmd Xl 

= (-±=\ y , /((k-p)/ e *)e^- k ) x 5 p , Xl dpdx i y'F(k,m)e jmx dm = ^ kx a kx , (96) 

2) fiF 



l y=\ J m • V p /((k - P )/e*)e i (P+ ac - k ) x F(k, ae)o P)Xl dpdmdxi 

^=) I /((k - p)/e*)im • xe i(p+m - k) ' x F(k, m)a P)Xl dpdmdx 1 

-j=\ j /((k - p)/s*)e iip+m - k) -*F(k, m)m • V p a p , Xl dpdx 1( im 
x • V x a; kx a kx + zV x u; kx • V k a kx , (97) 



34 



3) f F] 



(-j=\ J /((k - p)/ £ *y(P+ m - k )- x (p - k + m/2) • V k F(k, m)ap, xl dpdmd Xl 

= -i (-!-] J /((k - p)/s*)e iip - k) -*i(p - k)a p , Xl dpdxi • y e* mx V k F(k, m)dm 

+ Pf=\ J /((k - p)/ £ *)e i (P- k > x a p , xl dpdx 1 (-i) y i/2e imx m • V k F(k, m)dm 

= -iVk^kx • V x a kx - V 2 ( v k • V x )a; kx a kx , (98) 



4) fiF] 



j m • V p /((k - p)/e*) e i(p+m ~ k) ' x (P - k) • V k F(k, m)fip )X1 dpdmcki 

= / /((k - p)/e*)ie i(p - k) - x a p , Xl x • (p - k)dpdxi I -e imx m • V k F(k, m)dm 

+ / /((k - p)/e*)e i ( p - k >*a P) x 1 dpdx 1 I -e imx m • V k F(k, m)dm 

+ y /((k - p)/e*)e i ( p - k )- x V p S p , xl • (p - k)dpd Xl y -e imx m • V k F(k, m)dm 

= «(V k • V x )w fc:E x • V x a kx + i(V k • V x )a; kx a kx + (V k • V x )w kx (V k • V x )a kx . (99) 



B Calculation of £ 



In order to calculate £, we solve Eq. (78). Notice that cosh 2^ > 1 since [x ev > v. Using the definition of the 
cosh function, we can write 

e 2i + e -2i = 2 »ev 



and denoting £ = e 2 ^, one obtains a quadratic equation for £ with two solutions 

fJ-ev ± M 



£ 
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After we take into account Eq. (79), we obtain the following expression for £ 



e = -in- Mew 



4 fi ev + zv 



(100) 



C Bogolyubov transformation of the 0(e) part 



r(l) 
7,e' 



Here, we show how Bogolyubov transformation works on H f ^. We consider the terms of Eq. (71) starting 
with 



&(x • V x ^)a* H — (x • V x i^)(aa_ + a*a*_ 



dkdx 



/ («6 + t;6l)(x • V x/ u)(u&* + t;&_)dk<ix + - / (x ■ V^v)(ub + vb*_)(ub- + vb*)dkdx 
/ 6 [(u 2 + t> 2 )(x • V x ^ e „) + (u 2 - u 2 )(x • V x/ u od ) + 2nw(x • V x z/)] 6* dkdx - 
- f(uv(x ■ V x m) + ^(m 2 + v 2 )(x • V x i^))(66_ + b*b*_)dkdx 



+ c.c. 



m 2 - I — L v 2 



-6(x • V x c^)&* + ^ x • V x J 1 - — (**- + &*&-) 



dkdx. 



Here we used the following equalities 



u 2 + v 2 = cosh(2£), 2uv = sinh(2£), u 2 - v 2 = 1, 
sinh(2£) = d v u, cosh(2£) = d^u, 1 = <9 Mod w, 

2 / 2~ 

sinh(2^)VxM«, + cosh(2£)V x i/ = -^V X W 1 - -^. 

" V ^ 

Next, we consider the terms of Eq. (71) with the Poisson bracket 



id{fi,d*}dkdx-\ — / i{u, d-}dkdx = - \ ( a{fi, a*} + d{u, d_} )dkdx + c.c. 



/ (ub + vb*_)({fJb, ub* + t> &_} + {z/, u6_ + vb*}j dkdx + c.c. 



Note that 



(101) 



(102) 



(103) 
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• bb* terms give zero because their coefficients are purely imaginary and we add c.c. values in the end, 

• 6*V6 terms can be obtained as c.c. of bSJb* terms. 

Here, V denotes a gradient either with respect to x or with respect to k. Now, let us consider (103) term 
by term. 

1. bVb* and b*Vb: 



u 2 b{fi, b*} + v 2 b*_{n, &_} + vub*_{v, b-} +uvb{v, b*} 

v ' 

change k — ► — k 
I (b{u, b*} - b*{uj, b}\ dxdk = i f b{cu, 6*}dxdk. 



dxdk + c.c. 



Here, we used the fact that J b{u), b*}dkdx = — j b*{cu, b}dkdx 



(104) 



2. bb- and b*b* 



ubb-{/i, v} + vb*_b*{[i, u} + ubb-fa, u} + vb*_b*{v, v} 

V v ' 

give because v is even 



dkdx + c.c. 



l -j [(bb--b*b*_)(u{n,v}-v{v,u}) 

(bb_-b*b*_){nod,Z}]dkdx. 



dkdx = — 
2 



(bb- - b*b*_) {fj, ev + n odl £} 

V v ' 

only fj, oc i survives. 



dkdx 



(105) 



Here, we have used Eq. (100). 



3. bVb- and b*Vb*_: 
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uvb*_{fi, b*} + uvb{^, b-} + u b{u, &_} + v b*_{v, b*} 



tikdx + c.c. 



6sinh(2^){// eu ,6_} + 6cosh(2^){i/,6_} 

dkdx + c.c. 



dkdx + c.c. 



/J'ev 



H'ev 



We used Eq. (102) here. 



Finally, the rest of 0(e) terms are 



— / z>aa_dk(ix + c.c. = — / vbb-dkdx + c.c. 



(106) 



(107) 



Combining Eqs. (101), (104), (105), (106), and (107), we obtain the O(e) part of the Hamiltonian given by 
Eq. (82). 



D Canonicity conditions for near-identity transformation 

Here, we obtain the canonicity conditions for the coefficients of the near- identity transformation (83). For 
the canonicity up to O(e) order, we use the equation of motion in the Hamiltonian form. 



ib u 



6H_ 



:io8) 



In this Appendix for simplicity of notation, we skip writing x in the subscript of the dynamical variables. 
Using Eq. (83), we obtain 

*(c k + a k cl k + Ac{7k, c! k }) = — . (109) 

Since we are neglecting terms of the order higher than O(e), we can use the following approximation 



. SH 

{ 

Sc-u 



zw_ k c_ k . 



(110) 
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Combining Eqs. (109) and (110), we obtain 



Yy = J-* ~ «kw-kcl k - /3k{7k, w_ k cl k }. (Ill) 



Further, we have the chain rule in the form 



5 2L- f( 5H6 K | sh sb-* \ 



112) 



Using Eq. (83), we find 



K ~ k ' 

-^ = a- q ^-/9- q {7-q,5k}q' 

where the subscript of the Poisson bracket indicates the differentiation with respect to q. Therefore, Eq. (112) 
becomes 

YT = TIT + "-k^-kcl k + {7-k, /5-k^-kcI k }k- (113) 

5c k 5b k 

Combining Eqs. (Ill) and (113), we find 

= -2u_ k c\a od - /3 k {7k, w_ k cl k } + {7-k, /?-kW-k<i k }. (114) 

Finally, we obtain the canonicity conditions given in Eq. (84). 



E Near-identity transformation 

Lets us first apply the near-identity transformation to the 0(1) part of the Hamiltonian that is given by 
Eq. (80) 

/ cob*b dkcbc = / ujc*c + (ujo*cc- + oj/3*c{7*,c_} + c.c.)c?kdx + h.o.t. 
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To apply this transformation to the 0(e) part we just need to substitute b with c in Eq. (81). The non- 
diagonal terms cancel in the Hamiltonian if 

f (a + uja*)cc- + fu;/rc{7*,c_} + ^c{<^, c„}\ dkdx = 0. (115) 

Let us choose 

7 = 7 = ip = . (Ho) 

Then we can rewrite Eq. (115) as 

/ (<r + uja*)cc- + ( (u; ev + L0 o d)l3* + -r^- J c{(p, c_}dkcbc = 0. 

Integrating by parts one can show that 

/ uj od f3*c{ip,c^}dkdx = - {LL> od /3*,ip}cc-dkdx. 

Therefore, the diagonalizing condition becomes 

/ ( a + ua* + -{ui od j3*, if} J cc- + f u ev /3* + — ^ J c{cp, c_}dkcbc = 

From Eq. (117), the condition on (3 immediately follows 

/? = ^k (ii7) 

In order to obtain the condition on a, we expand the second term in the integral 

ua* = u od a* od + ^evOtlv + u od a* ev + uj ev a* od . (118) 

Integral over the last two terms vanishes because these functions are odd. Therefore, we consider only the 
other two terms. Next, we insert this expansion into Eq. (117) 

/ f <? + w od a* od + u ev a*e V + 2 U od{P*, V?} + -^P*Wod, </?} ) cc - + ( ^evP* + -^ J W, c_}cdkdx = 0. 
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Then, we obtain the following diagonalizing conditions on a 



OLev = ~ < ^od, > , (119) 

a od = -UP,—)- (120) 



Let us prove that a oc i = 0. Substituting Eq. (117) into Eq. (120), we obtain 

i ( n\ v u ev \ 
4 [uoj ev fi ev ) 

Expanding the Poisson bracket we find the following identity 

H 2 ev u ev \ u{fX ev ,UJ ev } + fJ,{u ev , v] + u> ev {v, fl ev } 



2 



:i2i) 



According to the definition, oj\ v = fi^ v — v 2 . Differentiation of both sides of the last equality with respect to 
x or k yields 



Vu ev = —SJ[i ev - —Vv. (122) 



iVev bJev 



We use Eq. (122) to show that 



{fi ev ,UJ ev } — {fj, ev ,v}, 

{u ev ,i/} = -^-{{i ev ,v}. 

Uev 



Plugging Eq. (123) into Eq. (121) proves that a d = 0. 
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